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I. INTRODUCTION 



Confining gauge theories generally become deconfined at high temperatures. The decon- 
fining phase transition in QCD is the subject of intensive experimental and lattice studies. 
While some general features of this phase transition can be extracted from symmetry argu- 
ments the details cannot be studied analytically due to strong coupling effects resulting 
from the non-perturbative nature of the phase transition. Therefore it is instructive to inves- 
tigate theories where deconfinement happens in the weak-coupling regime and hence can be 
studied quantitatively. One example of such theory, the (2+l)-dimensional SU(2) Georgi- 
Glashow model, is considered in this paper. The model consists of a SU(2) gauge theory 
coupled to a scalar (Higgs) field in the adjoint representation with a non-zero VEV v. By 
studying this model one may hope to find certain features of the deconfining phase transi- 
tions which this model shares with QCD. At zero temperature the (2+l)-dimensional SU(N) 
Georgi-Glashow model could be related through dimensional reduction to high temperature 
SU(N) gauge theory in four dimensions @, ||. 

It has been shown by Polyakov |3| that 't Hooft-Polyakov monopoles || [which are 
instantons in (2+l)d] give rise to the area law for Wilson loops, i.e., confinement. From 
the symmetry properties of the Polyakov loops |[| it follows that the Georgi-Glashow model 
must possess a sharp deconfining phase transition at some temperature. The first analytical 
attempt to find the temperature of the phase transition in the limit of the small gauge 
coupling g < t) was made by Agasyan and Zarembo 0. Following [| they argued that at 
low temperatures (T -C M w ) the theory is equivalent to a 2-dimensional gas of monopoles 
interacting with each other via (almost massless) Coulomb photons. (Remember that the 
original SU(2) symmetry is broken to U(l).) The phase transition in the monopole gas would 
occur through the binding of monopoles and anti-monopoles into pairs similar to Berezinsky- 
Kosterlitz-Thouless (BKT) phase transition |7j. From this picture they found the critical 
temperature to be T c = g 2 /(2n) ||, where g is the dimensionful gauge coupling. A criticism 
of this treatment was subsequently raised by Dunne, Kogan, Kovner, and Tekin ||. They 
argued that the correct description of the low temperature limit of the Georgi-Glashow 
model should include W bosons interacting with each other and with the monopoles by 
exchanging photons. The model thus reduced to a two-dimensional gas of point charges 
(Ws) and monopoles. The presence of thermally excited W bosons decreased the critical 
temperature by a factor of two yielding T c = g 2 / (4tt) ||. The authors of || also showed that 
the W's are required for the phase transition to be of the Z2 universality class, in agreement 
with symmetry arguments 

In this paper we re-examine the calculation of the critical temperature. We shall show 
that it actually depends on both the gauge coupling g and the ratio between the Higgs mass 
run and the gauge boson mass mw- Specifically, our result reads 

T. = f^, (1) 
where e = e{mn j mw) is a function of the ratio of Higgs and W masses M H and M w defining 



the monopole action 0, ITC 



9 \m w , 

When mil mw, which is the limit where the Bogomolny-Prasad-Sommerfield (BPS) 
bound on the monopole action is saturated, e(0) = 1 and T c in Eq. ([]]) coincides with the 
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value found in Ref. [§]. Outside this limit, our result disagrees with both Refs. || and 0. 

Our derivation of Eq. (|l|) is based on what we believe to be the correct interpretation of 
the renormalization group approach initiated in Ref. ||. To make the paper self-contained, 
we will provide alternative derivations for many facts already known in Refs. || |8| . Readers 
already familiar with Refs. |], || can jump directly to Sec. [V], which is the core of our paper 
where Eq. is derived. 

The paper is organized as follows. In Sec. |TJ we formulate the Georgi-Glashow model, 
mostly in order to introduce notations. In Sec. [ITI| we rewrite the partition function in term 
of a 2d Coulomb gas of electric and magnetic charges. In Sec. [TV] we show that the Coulomb 
gas, in its turn, can be recast into a field theory. We give two equivalent representations of 
this theory in terms of bosonic and fermionic variables. Our bosonic field theory is different 
from the one used in ||. Sec. [V] is devoted to the renormalization group analysis of the 
field theory where we determine the value of the critical temperature. Comparison with 
earlier works and a simple physical argument rederiving Eq. ([]]) are given in Sec. [VT|. In 
the Appendix we review the free massless scalar in two dimensions and show how a bosonic 
theory for the gas of electric and magnetic charges is constructed. 



II. THE MODEL 

In this paper we shall work mostly in Euclidean space. The Georgi-Glashow model is a 
SU(2) gauge theory where the gauge symmetry is partially broken by a triplet scalar, 

l = \{d^Y + K\<P\ 2 -v 2 ) 2 + \fI u (3) 

In the unitary gauge where <p l = <fi 2 = 0, the photon field A 3 is perturbatively massless, 
while W = A 1 ± iA 2 are massive gauge bosons with the mass mw = gv and are charged 
with respect to the photon field. The theory contains two dimensionless parameters: g/v, 
which is the expansion parameter of the perturbation theory, and m^/m^ ~ vX/v. We 
shall assume g/v <C 1 and mn/mw parametrically of order one, that is my/ ~ m# ^> g 2 . 

The model contains a classical solution — the 't Hooft-Polyakov "monopole", which 
is really an instanton in (2+l)d. (We shall use the terms "monopole" and "instanton" 
interchangeably.) At large distances, in the unitary gauge, the U(l) field of the instanton is 

where e is the magnetic charge of the monopole, 

47T 

e = 7 . (5, 

In our case ge is twice larger than the Dirac quantum 2n. The monopole action is given by 
Eq. (HD and is large at weak coupling g/v <C 1; correspondingly the density of instantons is 
exponentially small, ~ e _5 °. Due to the contributions of these instantons to the partition 
function, photon acquires a mass m 7 which is exponentially suppressed, m 7 ~ e~ s °^ 2 . 

In the Euclidean formalism, one puts the field theory at a finite temperature T by making 
the Euclidean time direction periodic with the period (3 = 1/T. It turns out that the 
interesting physics of deconfinement occurs at temperature of order g 2 , which is much smaller 
than W and Higgs masses but is much larger than the photon mass m 7 . We shall concentrate 
on the behavior of the theory at such temperatures. 
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III. COULOMB GAS REPRESENTATION 



Our goal is to rewrite the partition function of the Georgi-Glashow model at finite tem- 
perature T such that m 7 « T « m-w ~ m# in terms of an effective two-dimensional 
theory. Extending Ref. Q, our first step is to rewrite the partition function in terms of a 
two-dimensional Coulomb gas of monopoles and W bosons, 

oo 1 . M N 

z -Ej^EU ^ II (6) 

M,AT=0 ' ' ei,ffj t=l i=l 

where ^(xj, e^; jjj, gj) is the effective action of a system of M monopoles located at the (2d) 
positions X{ with magnetic charges e*, and A 7 " IV-bosons located at positions ?/} with electric 
charges gj . The sum over e« , ^ in Eq. (|j) goes over different signs of the W-boson and 
monopole charges. The effective action is the sum of three contributions, 

S(Xi, &i\ 2/jr , gj) = S mon (Xi, Si) + SwiVji 9j) + S'mt( x ii e u Vji 9j) (7) 

where S mon is the part of the action coming solely from the monopoles, Sw is coming solely 
from the Ws, and S- m t represents the interaction between the monopoles and the Ws. 

The explicit form of each term in Eq. (|7|) can be found from physically intuitive argu- 
ments. For S^on, we notice that the monopole gas is effectively two-dimensional because our 
assumption about the temperature T ~ g 2 ^> m 7 ~ e ~ s o/^ ma k es the size of the compact- 
ified time direction (1/T) much smaller than the average distance between the monopoles 
(~ e 5 °/ 2 ). Thus, the monopoles interact via the 2d Coulomb potential p, and 

e 2 T 

S m0 n(xi,ei) = MS - — y^ejej ln(T|fj - xj\) , ^ei = 0. (8) 

ij i 

Here = ±1 are the signs of the monopole charges, which are measured in the units of e 
defined in Eq. (0). The action in Eq. (^) integrated over all monopoles spatial positions is 
finite only when the total magnetic charge vanishes; otherwise it is logarithmically divergent. 
The reason why T appears in the argument of the logarithms in Eq. (^) is that T" 1 is the 
smallest distance where the inter-monopole potential is logarithmic. 

Turning to W bosons, we notice that since T <C Mw the Ws are almost static. As they 
carry U(l) electric charges they also form a Coulomb gas. We can visualize the worldlines 
of the W bosons as small circles wrapping around the Euclidean time direction; at distances 
large compared to (3 = 1/T, the ensemble of such circles are indistinguishable from a gas 
of point-like objects. The action Sw is simply the exponent in the Boltzmann suppression 
factor, i.e., Ew/T where E\y is the energy of the W configuration. Therefore for almost 
static Ws 

2 

Sw(Vj, 9j) = N^- - — 9i9j H T \Vi -Vj\), 9i = ±1 > ^9i = • (9) 

ij i 

The first term in the right hand side comes from the static energy, and the second term from 
the interaction between the Ws. Although the contribution of W bosons to the free energy is 
exponentially suppressed by e~ mw ^ T , they can still destroy confinement, since confinement in 
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this model is due to instantons and the density of instantons is also exponentially suppressed 
(0, also see below). 

The origin of the interaction term S- m t is less trivial: it is due to the complex phases 
acquired by the W bosons along their trajectories (which are closed time loops) in the U(l) 
field created by the monopoles. This purely imaginary contribution to the effective action 
of the Coulomb gas has been explicitly computed in ref. to be 

N M 

-Sint = 2i 2J ^ eigj9(xi - yj) (10) 
t=i j=i 

where 8(xi — y*j) is the angle between the vector connecting the monopole at Xi and the 
W boson at yj and a chosen spatial direction (say, the x 2 axis). The sum in Eq. JIHD does 
not depend on the choice of the reference direction if the gas is electrically or magnetically 
neutral, or both. Thus it is always well-defined when S mon and Sw are finite. 

Equation flID"D can be derived in a simple, almost geometrical, way. Let us consider a 
simple case of one monopole and and a pair of W + and W~ bosons. First we "unwind" the 
time direction, so that we have W bosons move along finite time intervals in the field of 
an infinite number of monopoles, positioned periodically along the time axis. It is obvious 
from Fig. |l] that the phase obtained by a W boson in this case is the same as when it moves 
along an infinite temporal trajectory in the field of a single monopole. 

M 



P 



M 
M 
M 
M 



W 



= MO 



W 



FIG. 1: The phase obtained by a IF in the field of periodic instantons (see text). 



Now the total phase obtained by a pair of W + and W~ bosons is proportional to the 
magnetic flux passing through the two-dimensional surface stretched between the worldlines 
of the two bosons (see Fig. |2|). If we surround the monopole by a sphere, then the flux 
is proportional to the area between the two meridians as shown in Fig. Projected to 
the equator plane, one sees that flux is indeed proportional to the angle between the line 
connecting the monopoles to the W + and W~ bosons, in agreement with Eq. fllcf) . The 
coefficient 2 in Eq. flID"D is due to the fact that the U(l) charge of W bosons is twice the 
elementary value. 

Gathering all terms, we arrive to the Coulomb-gas representation of the partition function 
of the Georgi-Glashow models at temperatures large compared to m 7 , but small compared 
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w~ 




FIG. 2: Geometrical interpretation of the phase obtained by a pair of W + and W boson in the 
field of a monopole. 



to mw, 



z =n%% e e m^n^-p 

MAT ' ' e,=±l ffj=±l J i=l 



N 



3=1 
M N 



2vr 



2ttT 



i=l 3=1 



11^ 



Here Co an d Co are the instanton and W fugacities, respectively. They have dimension 2 and 
are both exponentially suppressed, 



">2 — So 



Co ~ TV 
Co ~ T 2 e~ mw/T . 



(12a) 
(12b) 



(The pre-exponents will not be needed for our computation of T c .) We put the subscript "0" 
in the notations because these fugacities serve as the initial conditions for the renormalization 
group in Sec. [V]. The sums over e, and are restricted to configurations with zero total 
electric and magnetic charges. 



IV. EFFECTIVE FIELD THEORIES 

In this section we (re) derive effective bosonic and fermionic theories describing a 2d gas 
of electric and magnetic charges. Our bosonic Lagrangian is different from the one obtained 
by the authors of || . 
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A. Bosonic theory 



By using the result of Appendix A, one can write each term in the sum ( [LTD as a multipoint 
Green's function, 



M\ JV' ^ ^ 

MN ' ' ei=±l ffj=±l 



where we have introduced 



N 



\i d 2 x t n ^ ( t n e*^*^ n e ^ 



i=l 



i=i 



/5 = eVT, /3 







T 



(13) 



(14) 



and the average is taken in the free massless bosonic field theory with the ultraviolet cutoff 
at the scale T. $ is a (scalar) bosonic field and is the dual field. 

One can now perform the summation over e, and Qj and over M and N in Eq. fll~3|). The 
result reads 

d 2 x ( 2Co cos /3$ + 2Co cos /30 N 



Texp 



(15) 



The right hand side, according to the well-known formula of perturbation theory, is equal 
to 

Tr e~ HT 

Z= _ 6 (16) 



Tr e~ H ° T 



where T is the total Euclidean time and 



H 



dx 
dx 



l -{d x §) 2 + 1 -{8 X Q) 2 + 2C cos/3$ + 2C cos/36 



(17) 
(18) 



Thus, our problem is reduced to finding the ground state energy of the Hamiltonian H from 
Eq. ([TJD, defined with an ultraviolet cutoff at T. 

One can construct a (Euclidean) Lagrangian corresponding to the Hamiltonian of Eq. (|18[) 
by requiring that the free part of the Lagrangian gives the same correlation functions as in 
Eqs. ([AT]) and (|AJ). The result is 



-(a,$) 2 + -{d x ef - id x <& d T Q + 2Co cos /?$ + 2C cos (30, 



(19) 



where $ and O should now be understood as two independent fields. Note that the La- 
grangian in Eq. (|19"D is different from the one derived in || . The main difference is that our 
Lagrangian ( p~9|) is written for two independent scalar fields and the Lagrangian used in || 
employs a scalar field x an d its dual field x satisfying the duality condition 



duX 



-ze 



fJtU 



(20) 



However, the duality condition (|20|) implies d 2 x = d 2 X = 0, that is the dual field in the sense 
of Eq. ( pOD exists only for harmonic functions. It is therefore difficult to define Feynman 
functional integrals over x f° r t ne partition function employing the Lagrangian which uses 
both fields x an d X since the integration would have to be restricted to harmonic functions 
only. 
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Defining a vector field 



C M = dfl - ie^d u e (21) 
the Lagrangian of Eq. ([H|) could be rewritten in a more Lorentz (rotational)-invariant form 

£b = ^(^e) 2 + ^(n M C At ) 2 + 2Cocos/?$ + 2C cos/3e, (22) 

where n M is a unit vector in a randomly chosen direction in the (r, x) plane. Putting 
= (0, 1) one would obtain the Lagrangian in Eq. ([TDD - 

In the limit of no monopoles we put Co = in Eq. fl2~2"|). Integrating out the $ field yields 
us the Lagrangian of the sine-Gordon theory for the 2d gas of W bosons 

£Jg = ^e) 2 + 2C cos^e. (23) 

Similarly putting (q = in Eq. ( p2|) would eliminate W bosons from the theory and would 
give us a sine-Gordon theory for monopole gas, similar to || 

tic = i(^$) 2 + 2Cocos/3$. (24) 

The exact solution for the S'-matrix of the sine-Gordon theory is known [11| and the exact 
expression for the soliton mass has been found |f[2]| , allowing for quantitative study of the 
limits shown in Eqs. (^) and ( 24|) along the lines of |J. An exact solution of the full 
theory given by Eq. ( p2j) has not been found. This will prevent us from making quantitative 
predictions about the behavior of the string tension around T c but will not interfere with 
our determination of T c itself. 

Polyakov loops play the role of the order parameter of the deconfining phase transition 
To compute the correlation function of two Polyakov loops, one inserts two fundamental 
charges into the Coulomb gas ([□]). This means that Polyakov loop P(x) is mapped to the 
operator e l/3e ^ 2 of the effective theory, so that || ^[ 

(P(f)pt(o)) ~ ^exp (^e(f)) exp (-^ (°)) ) ( 25 ) 
where the average in the right hand side is taken in the theory described by Eq. (18). 



B. Fermionic theory 



It is easy to apply the bosonization rules (T3| to map the bosonic theory (|P8|) onto a 
fermionic theory. We noticed above that without the Co cos $® term the Lagrangian of 
Eq. fl22|) reduces to a sine-Gordon theory given by Eq. (|24]) . As was shown in the 
sine-Gordon theory is equivalent to_ the massive Thirring model with Dirac mass term. 
Mandelstam's rules [|X3j] map Co cos fi® onto a Majorana mass term. Thus, the fermionic 



effective theory is given by the Lagrangian 



— G — — Tfl — — 

C f = ip^d^ + -(VW) 2 + + -(^ + W) (26) 



S 



where ip c — Cip T is charge-conjugated to ip. In the chiral basis 

*=(£)• *=(-5)- (27) 

The four-fermion coupling G is related to the parameters of the bosonic theory by 

1 + ° - 471 - P (28) 

1 + 7"^-4^' (28) 



while the fermion bare masses are 



Co ~ Co / ori N 

m~— , m~— . (29) 



One can check by an explicit calculation that the Lagrangian ( Pq) gives the partition function 
ofEq. ©• 

Defining real Majorana fermion fields M 

P=^f, Q = ^f (30) 
we can rewrite the Lagrangian of Eq. (|2~6D as 

C f = P-y^P + Q-y^Q - 2GPPQQ + (m + m)PP + {m-m)QQ. (31) 



The physical meaning of the Lagrangian in Eq. (plD will be clarified in the next Section. 

The bosonic theory (fL9|), and, therefore, the partition function, is symmetric under duality 
transformation Co Co an d (3 «-> j3, or, equivalently under $ <-> O. The duality transforma- 
tion $ <-> corresponds to ip% *-> ipi, ip2 ipt- This transformation interchanges the Dirac 
and Majorana mass terms, and the U(l) and U(1)a currents j M = ipj^ip and j 5tl = ^j^^ip. 
Naively under this transformation the interaction term changes sign, G <-» — G, which is 
not consistent with (3 <-> (3 in Eq. (p8|) . This is because the symmetry between the Dirac 
and Majorana mass terms in Eq. (|26|) is broken by the implicit procedure of regularization: 
when computing the loop graphs in the theory (|26|) one should use a regularization scheme 
which preserves the conservation of when rh = 0, e.g., the Pauli-Villars scheme with a 
fermion with large Dirac mass Jl4]]. Using a different regularization scheme, for instance 



the Pauli-Villars scheme with a heavy Majorana fermion, would result in a relation between 
fermionic and bosonic couplings different from Eq. (ESI) fT4||. 



V. RENORMALIZATION GROUP, MASS GAP AND CRITICAL TEMPERA- 
TURE 

A. Review of RG in sine-Gordon theory 

At any temperature, the correlation length of the Georgi-Glashow model is equal to 
the inverse mass gap in the effective theories (p2|) and (^6|). The fact that there are two 
perturbations on top of the free bosonic theory in Eq. ( f22|) is crucial for the existence of a 
critical point with infinite correlation length (zero mass gap). To understand the mechanism 
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of the disappearance of the mass gap, let us recall how the gap emerges when there is only 
one perturbation. 

We suppose, for a moment, that Co = 0. Then our theory is reduced to the sine-Gordon 
(SG) model of Eq. (^4|). The full spectrum of the SG theory is known HIT], |I2~|, but for 



24]). The full spectrum of the SG theory is known fTI , 
our purposes we only need a method for a parametric estimate of the gap. The Wilson 
renormalization group (RG) can be employed to this end. In this procedure we decrease 
the ultraviolet cutoff A by successively integrating out the modes with momenta above the 
cutoff, but below the value of the cutoff at the previous RG step. Since the fugacity Co will 
run with A we will denote the renormalized fugacity by the function ((A). 

When ( is small, the perturbation is small and the RG equation for ( is set solely by the 
conformal dimension of the operator cos /3$ jl4j> [H| 



d\ A 2 1 j A 2 ' 



A 



i T 
ln A 



A = — 



47T 



which has an obvious solution 



_C_ 

A 2 



J^2 



2-A 



47T 



(32) 



(33) 



We have written the RG equation at the leading order in the dimensionless coupling C,/A 2 
which controls the conformal perturbation theory. If (3 2 > 8ir the perturbation is irrelevant 
since C/A 2 is smaller than ( /T 2 for all A < T. In the case of j3 2 < 87r where the perturbation 
is relevant, there exists a scale where C/A 2 ~ 1 and the perturbation is no longer small. Our 
RG equation (^) is not valid beyond this point. The theory becomes strongly coupled at 
this scale and one can expect the emergence of a mass gap. Thus the mass gap is of the 
order 

f \ V(2-A) 

\ / (||) ^ Te -5o/(2-A) (34) 



11, 12 



which agrees parametrically with the result of 
not allow the reproduction of the pre-exponents in Ref. 
the exact expression for the ^-matrix of the theory. 



The nature of the RG method does 
for which one actually needs 
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B. Renormalization group and criterion for T c 



Now we consider our case, where both perturbations are present [see Eqs. ( |18l) and (gg)]. 
We now have two RG equations 



with solutions 



AH 

d\ a 2 
d\ a 2 



A 2 



(2-A) 
(2-A) 



T 2 1 A 



2-A 



c 

A 2 ' 


A 


_ e 2 T _ 

47T 


4vrT 
9 2 


c 

A 2 ' 


A 


9 2 

4nT ' 






c 


Co (T\ 


2-A 




A 2 


T 2 \A J 





(35a) 
(35b) 

(36) 



Noting that A A = 1, one can immediately distinguish three temperature regimes. 
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(i) T < g 2 /(8%). In this case A < ~, A > 2. The operator cos /3Q is irrelevant and can 
be ignored. The infrared properties of the theory is the same as of the sine-Gordon 
theory fl2~4|) and hence the theory has a mass gap || [12]. 



(ii) T > g 2 /(2n). In this case A > 2, A < |, the operator cos/3$ is irrelevant and the 
theory is in the massive phase of the SG theory of Eq. (p3[). 

(hi) g 2 /(87r) <T < g 2 /(2ir). In this case both perturbations are relevant. 

The critical point should be in the regime (iii) which we now consider. If we follow the 
Wilson RG procedure, then at first, when A ~ T, both perturbations are small. As one 
decreases A both perturbations become more and more important. However, if at some 
scale one perturbation is of order one while the other is still small, then one can_ conclude 
that the theory has a mass gap. Indeed, suppose that at some A, (/A 2 ~ 1 while (/A 2 <C 1. 
In this case the theory is a SG model with a very small perturbation ~ (. Since the SG is 
fully gapped, the perturbation theory over C/A 2 is non-degenerate and does not destroy the 
existence of the gap. The same happens when (/A 2 becomes ~ 1 when (/A 2 is still small. 
Therefore, the theory can be critical only when (/A 2 and (/A 2 are of order one at the same 
energy scale. 




FIG. 3: The running of fugacities at T = T c . 



In Fig. |3] we illustrate what happens at T c . At the scale T = T c we chose the monopole 
fugacity ( to be much smaller than the W fugacity (. (In principle it can be the other way 
around as well.) However, ( runs faster than £, so that at some scale \i the two fugacities 
become equal and are of the order /i 2 . In other words, the perturbative expansion over the 
two fugacities breaks down at the same scale \i where C/fi 2 ~ C/^ 2 ~ 1- We can guess that 
when this condition is satisfied the gap disappears. 

If our guess is correct, then T c can be found be requiring the existence of a scale /i where 

2 1. (37) 



Using Eq. ([TJ) in Eq. (|37]) yields 



e -*(I) ~^g) -I- (38) 
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Solving this equation, taking into account Eqs. (|2|) and (14]), we then find our final result 
for T c , 

T - g2 e + 2 
c 4vr 2e + 1 ' 

and an estimate for the non-perturbative scale /x at the critical temperature 

^-Texpf —J, T^T C . (39) 

We notice, however, that the RG procedure above does not prove that the temperature 
found above is the critical temperature: it merely says that if the theory becomes critical at 
some temperature, then this temperature should be given by Eq. ([!]). 

To argue that the mass gap vanishes at T c the fermionic Lagrangian of Eq. ([31]) becomes 
useful. This can be seen most clearly in the limit of the light Higgs, m# <C m^ ||. 
(We still have ran 3> g 2 , or, equivalently, mn/mw g/v, so that one loop corrections 
to the monopole action @ would be finite and Eq. @ would still be valid |L6|.) Then 



T c = £ 2 /(4vr) with (3 2 = (3 2 = 4vr, which leads to G = [see Eq. fl28|)]. The Lagrangian (|3T]) 
becomes that of a theory with two free Majorana fermions with masses mp = m + m and 
rriQ = \m — m\. Using Eq. fl29|) with the critical temperature condition Eq. ([37]) we may 
argue that at T c one may get m = rh. When m = rh we obtain mp = 2m and mQ = 0. One 
of the Majorana fermions is massless, and the other one is massive. The existence of the 
massless Majorana fermion near the phase transition is in agreement with the Ising nature 
of the phase transition 0. Notice that, strictly speaking, the RG equations can give us 
mp~/i, but cannot show that mg = 0. In the bosonic language, there is no mass gap when 

/3 2 = /3 2 = 47randCo = Co & 

When the ratio mn/mw is large, the fermion coupling G is also large. While we do 
not know the exact mass spectrum of the theory (j3T[) , one may guess that qualitatively 
the picture does not change. This means that near T c the theory contains two Majorana 
fermions: a very light one with some mass mg, and another one with mass mp~/i> m_. 
Because G / the two fermions interact with each other. However, at energy much below /i, 
the heavier fermion P decouples, leaving the dynamics to be completely determined by the 
lighter fermion Q. Since in 2d there is no relevant or marginal self-coupling of one Majorana 
fermion (e.g., {QQ) 2 = 0), the dynamics below \i is determined by one free Majorana fermion, 
as in the limit mH/mw — » 0. 



C. Correlation length, critical region, etc. 

It is instructive to find the behavior of the correlator of two Polyakov loops when T is 
very close to T c . This correlator has various behaviors at different length scales: 

(i) At very large distances \x\ 3> M _1 , where M is the mass gap in the theory, 

(P(x)P(0)) ~ e~ Mx . (40) 

(ii) In the intermediate range \x\ M _1 , where /1 is defined in Eq. ([39]), the 
correlator behaves that that of the order parameter in the critical Ising model, so that 

(P(x)P(0)) ~ |£|- 1/4 . (41) 
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(iii) At distances T 1 <C \x\ <C M 1 , the effective theory is the theory of one massless 
scalar field, and the behavior of the Polyakov loop correlator is determined by the 
conformal dimension of the operator e 4/30 / 2 , 

(P(f)P(0)> ~ |£|-(2e+l)/(2e+4)_ (42) 

It is interesting that this model has Ising-type behavior only at distances larger than /i -1 , 
which is a length scale exponentially large compared to the inverse temperature. At distances 
between T -1 and /i -1 the correlator of Polyakov loops has a power-law behavior, but with 
the critical exponent in Eq. ([|2J) which does not have an absolute fixed value and can vary 
continuously depending on m# /raw 

From the discussion above it also follows that the critical region, i.e., the range of temper- 
atures where the model exhibits the Ising-like behavior, is very narrow. For a temperature 
T to be in this region, the fugacities ( and ( should become large (~ A 2 ) almost at the same 
scale, where "almost" means a mismatch by a factor not much larger than one. This gives 

E^IA < 9 (43) 



T, r>y 
c V 

as the critical region. Deep inside this region, the correlation length exhibits the Ising 
behavior with the critical index v — 1: 

^\T-T c \-\ (44) 

Outside the critical region, the correlation length is set by the energy scale at which one 
of the perturbations cos/5<l> and cos/30 becomes large. Thus, 



( 1 / 2negv ^ 



1 / Angv 
— exp 



7^ ex P ( 77T7^ "2 ) ' T >T C . 



(45) 



One can see from Eq. (fl5]) that as the temperature increases, the correlation length first 
increases exponentially from m' 1 ~ e s o/ 2 a t small T to when T is near the critical 
region, and then decreases exponentially as T grows past the critical region. Inside the 
critical region itself, £ in Eq. ( fiST) has the same exponential factor as in Eq. (|39D but 
the prefactor of the expression for £ diverges at T c . Finally we note that in our model, the 
correlation length of Eq. (|45|) is exponentially large for all T mw ~ gv. 



VI. DISCUSSION AND CONCLUSION 

We have shown that the phase-transition temperature of the Georgi-Glashow model is 
given by Eq. (|l|). The critical temperature is of order g 2 and also depends on the ratio 
van / raw via the function e(raH / raw) defining the monopole action through Eq. @. In the 



Georgi-Glashow model e runs from e(0) = 1 to e(oo) w 1.787 0, [10| . If we formally take 
e — > limit in Eq. (H), which of course can never be physically realized in Georgi-Glashow 
model, we get T c = g 2 /(27i). This is the value of the critical temperature suggested in 



13 



Ref. ||. This is not surprising since when e — > the monopole fugacity becomes (o ~ 1 
while the W bosons' fugacity is still small £o <C 1, so that the W bosons are completely 
outnumbered by the monopoles and the problem reduces to physics of monopole-only gas 
considered in Ref. ||. In the opposite unphysical limit of e — > oo, the critical temperature in 
Eq. ([!]) becomes T c = g 2 /(87r). In this limit there is no monopoles and the phase transition 
temperature is simply that of the BKT phase transition for the 2d gas of W's. 

In the light Higgs limit g/v <C majm^ <C 1 corresponding to e = 1 our critical tempera- 
ture from Eq. ([I]) becomes T c = g 2 /(4n) in agreement with Ref. M. For other values of e, or, 
equivalently, for not extremely small values of the ratio tuh/tuw, our critical temperature 
is different from the one found in M . 

The methods used in our paper are almost identical to those of Ref. ||. The crucial 
difference is that in Ref. || the RG equation is run to A = ln(T/A) = oo, reaching what the 
authors of j5| interpret to be a fixed point at ( = ( = oo. In contrast, in our paper the RG 
equation is used only up to the point where one of the fugacities is of order one. Strictly 
speaking, the RG equations are derived at the leading order of perturbation theory in (/A 2 
and CM 2 , and thus cannot be used outside the regime (/A 2 , (/A 2 <C 1. 

There exists an intuitive argument leading to our result ([!]). According to this argument, 
the phase transition occurs when the densities of monopoles and W bosons are equal |§. 
Naively the instanton density is e~ So , but in 2d the instanton action acquires a logarith- 
mic contribution from distances between 1/T and the mean inter-instanton distance l mon . 
Therefore l mon can be obtained from the consistency condition equating the monopole den- 
sity obtained from geometric considerations ~ lAmon to the same density obtained from the 
modified monopole action. The condition reads 



1 exp 



I 2 



e 2 T 

-S - ——\n(Tl r , 

47T 



(46) 



Analogously, the mean distance between the W bosons lw is given by 

1 



— ~ T 2 exp 



/ 



m w g 2 , (rrl \ 
■ \n{Tl w ) 



(47) 



T AixT 

Solving Eqs. (EBf) and d4~7f) for l mon and lw, and requiring l m = lw, one then obtains Eq. (II]). 1 
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Almost the same argument was given in Ref. [pi, but the logarithmic terms in the exponents were missing. 
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APPENDIX A: FREE MASSLESS BOSON IN 2D 



We consider a free massless boson in two dimensions, 

S = j£x l -{d^f. (Al) 

Although Eq. (|A1|) is invariant under 0(2) rotations, we can choose an arbitrary direction to 
be the Euclidean time axis, and quantize the theory canonically. The canonical commutation 
relation reads 

[$(T,x),d T $(T,y)] = 6(x-y), (A2) 

and the Hamiltonian is 

h= 1 - jdx[{-d T $y + {d x <s>y\. (A3) 

We introduce a field 9 defined as follows, 



Q(r,x)=i dyd T $(T,y). (A4) 



Since r is the Euclidean time, is a Hermitian operator. By definition, d x Q = id T <&. 
One also finds d T Q = —id x <& by using the free field equation satisfied by $. Thus, d^Q = 
—ie^ v d v ^. We will say that O is dual to $. The commutation relation between $ and O 
can also be found from the definition of O 



[$(r,x), Q(T,y)] = i6(y-x). (A5) 



In other words, $ and O are mutually non-local. The Hamiltonian ( |A3[ ) can be written in 
a form symmetric under duality transformation $ <->• O, 

H= l -Jdx [(d x $) 2 + (d x Q) 2 ] . (A6) 

Since ( |A1|) is a free field theory, the Euclidean Green's functions can be easily computed, 

(T$(£)$(0)) = (T9(5)e(0)) = -— ln(m\x\) (A7) 
(T$(x)6(0)) = ± sgn(r) arccos = ^-9{x) (A8) 



where x = (r, x) and m is a small infrared regulator. In Eq. ( |A8|) the value of arccos is 
taken to be in (0, 7r). The function 9{x) is the angle between the vector x = (t,x) and the 
x axis, and runs between — 7r to n. The discontinuity of ( |A8| ) at r = 0, x < is due to the 
Euclidean time ordering, 

(T$(e, x)0(O, 0)) - (T$(-e, x)6(0, 0)) = ([$(0, x), 6(0, 0)]) = iB(-x), (A9) 

with here being the ^-function. The presence of a cut along the spatial axis means that 
the correlator (|A8|) is not invariant under 0(2) rotations in the Euclidean space-time. This 
is not surprising, since $ and 6 are not mutually local. 
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The following correlator can now be easily computed 



M n 
i=i j=i 



-((EiAJ'+CE^j) 2 )/^ 



exp 



M,N 



+ E^ ln (Ai^-%i)-^E#^ 



2tt 



(A10) 



where A is the ultraviolet cutoff. The m dependence tells us that the correlation function 
vanishes in the limit m — > unless Ylifti = Pj = (charge neutrality). The dependence 
of the correlator on the ultraviolet cutoff A is A _A , where 



47T 



(All) 



is the sum of the conformal dimensions of the operators in the correlator ( |A10| ). The 
correlator (|A10|) does not have branch cuts if all products (3i(3j are multiple of 2n; in this 
case, the result does not depend on the choice of the time axis. This is due to the fact 
that the operators in the correlator ( |A10|) are mutually local. For the correlator in Eq. ([13]) 
(3(3 = 47r and the partition function is indeed independent of the choice of the time axis. 
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